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Abstract: We study the poles of the retarded Green functions of a holographic supercon- 
ductor. The model shows a second order phase transition where a charged scalar operator 
condenses and a U{1) symmetry is spontaneously broken. The poles of the holographic 
Green functions are the quasinormal modes in an AdS black hole background. We study 
the spectrum of quasinormal frequencies in the broken phase, where we establish the ap- 
pearance of a massless or hydrodynamic mode at the critical temperature as expected for 
a second order phase transition. In the broken phase we find the pole representing second 
sound. We compute the speed of second sound and its attenuation length as function of the 
temperature. In addition we find a pseudo diffusion mode, whose frequencies are purely 
imaginary but with a non-zero gap at zero momentum. This gap goes to zero at the criti- 
cal temperature. As a technical side result we explain how to calculate holographic Green 
functions and their quasinormal modes for a set of operators that mix under the RG flow. 
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1. Introduction 

The AdS/CFT correspondence is a tool for studying field theories in the strong coupling 
regime [1-3]. The range of physical phenomena to which it can be applied is constantly 
increasing. One of the latest additions is the realization of spontaneous symmetry breaking 
and the appearance of a superfluid (often described as superconducting) phase at low 
temperature. A model with a charged scalar condensing in the background of a charged 
AdS black hole has first been introduced in [4]. Shortly afterwards it was realized that a 
charged scalar condenses as well when the black hole is neutral and it was shown explicitly 
that the DC conductivity is infinite in the broken phase [5]. 

By now there is a large variety of holographic models of superfluidity/superconductivity 
[6-21]. Specifically hydrodynamical behavior in these models has been addressed before 
in [22,23] where the speed of sound has been calculated from derivatives of thermodynamic 
quantities. The hydrodynamic poles of retarded Green functions have been studied in an 
analytical approximation for infinitesimal condensate in a p-wave model in [24] 

In this work we are interested in the hydrodynamics of the holographic superconductor 
introduced in [5]. In general, hydrodynamic behavior is connected either to the presence 
of a conserved charge, a spontaneously broken symmetry or a second order phase transi- 
tion. As we will see in our model all three possibilities are realized. In holographic models 
the hydrodynamic modes appear as quasinormal modes in the AdS black hole background 



- 1 - 



whose frequency vanishes in the zero momentum limit [25,26]. Quasinormal modes play an 
important role in the physics of black holes and a very good review of their properties in 
asymptotically fiat black holes is [27]. In the asymptotically fiat situation the quasinormal 
modes are defined as solutions of linearized wave equations with purely infalling boundary 
conditions on the horizon outgoing ones at the boundary. In asymptotically AdS spaces the 
situation is somewhat different. On the horizon one still imposes infalling boundary condi- 
tions but on the conformal boundary of anti de Sitter space several possibilities of choosing 
boundary conditions arise: Dirichlct or Neumann or Robin (mixed) ones. The holographic 
interpretation fixes this degeneracy of boundary conditions by defining the quasinormal 
modes as the poles of the holographic Green functions [28-30] . This of course implies to 
know the holographic Green functions, which are computed using the prescription given 
in [31]. In situations in which there are several fields whose linearized wave equations form 
a coupled system of differential equations possibly subject to a constraint due to a gauge 
symmetry the construction of the holographic Green functions is a bit more complicated. 
We solve this problem in full generality and show that the quasinormal modes defined are 
the zeroes of the determinant spanned by the values at the boundary of a maximal set of 
linearly independent solutions to the field equations. 

Having solved the problem of defining the holographic Green functions we concentrate 
on finding the lowest quasinormal modes and in particular the ones representing the hy- 
drodynamic behavior of the system. Hydrodynamic modes can be understood as massless 
modes in the sense that lim/t-^o ^(^) = 0- Such modes arise in the presence of a conserved 
charge. In this case a local charge distribution can not simply dissipate away but has to 
spread slowly over the medium according to a diffusion process. Other situations in which 
hydrodynamic modes appear are at a second order phase transition, characterized by the 
appearance of a new massless mode and spontaneous breaking of a global symmetry where 
a massless Goldstone boson appears. A discussion of hydrodynamics in systems with spon- 
taneous breaking of global symmetries can be found in [32] and in the relativistic context 
in [33]. 

We will consider the abelian gauge model of [5] without backreaction, i.e. assuming 
that the metric is a simple AdS black hole with flat horizon topology. In [5] it has been 
established that this model undergoes a second order phase transition towards forming a 
condensate of the charged scalar field thus spontaneously breaking the C/(l) gauge symme- 
try. The conductivity in the broken phase has a delta function peak at zero frequency and 
a gap typical of superconductors. We add here that although this model is referred to as 
holographic superconductor it is more proper to speak of a holographic charged superfluid 
as in [22] since the U{1) symmetry in the boundary field theory is global and there is no 
clear holographic description of how to add gauge fields. 

An important particularity of the model is the missing backreaction. Since the metric 
fluctuations are set to zero this means that effectively there is no energy momentum tensor 
in the field theory dual. In particular the generators of translations and rotations are 
missing in the operator algebra. This does of course not mean that the model does not have 
these symmetries, they are however not realized as inner automorphisms of the operator 
algebra (they are still outer automorphisms). Besides the space time symmetries being 



-2- 



realized as outer automorphisms there is a direct consequence of this in what concerns the 
hydrodynamics of the model: all hydrodynamic modes related to them are missing. There 
is no shear mode for the momentum diffusion and no sound mode for the energy transport. 
The hydrodynamic modes we find in the model are therefore only due to the presence of 
the C/(l) symmetry and its spontaneous breakdown. 

Hydrodynamics can be understood as an effective field theory defined by the continuity 
equations of the conserved currents and so-called constitutive relations which encode the 
dissipative behavior of the system. The constitutive relations tell us how fast a current is 
built up due to gradients in the charge density or due to external fields. The constitutive 
relations depend on transport coefficients such as viscosity or conductivity. Transport coef- 
ficients can be divided into absorptive and reactive ones depending on whether they are odd 
or even under time reversal [32]. A typical example for an absorptive transport coefficient 
is the diffusion constant, an example for a reactive one is the speed of sound. Reactive 
transport coefficients such as the speed of sound (the static susceptibility is another ex- 
ample) can often be computed from purely thermodynamic considerations. This has been 
done for the speed of second sound in this model in [22] and for a variant of fourth sound 
in [23]'''. As we have argued before, the hydrodynamic modes of energy and momentum 
transport, shear and sound modes arc missing. In the broken phase one expects however 
the appearance of a hydrodynamic mode with approximately linear dispersion relation for 
small momenta which represents the second sound present in super fluids. Indeed such a 
mode is bound to appear for each spontaneously broken continuous symmetry [32]. Our 
aim is to find the second sound mode directly as a pole in the holographic Green functions 
in the broken phase and to read off the speed of sound from its dispersion relation. We 
have found this mode numerically and our results for the speed of second sound agree (with 
numerical uncertainties) with the results in [22]. 

Below the critical temperature one expects actually only a part of the medium to be 
in the superfiuid state whereas another part stays in the normal fluid phase. The fluid 
is a two component fluid and naively one might expect that this is reflected in the pole 
structure of the Green functions as the presence of the diffusive pole of the normal fluid 
component. As we will see, the hydrodynamic character of this diffusive pole is lost however 
below the critical temperature. We find a pole with purely imaginary frequencies obeying 
a dispersion relation roughly of the form lj = —i^ — iDk^. The gap 7 goes to zero at the 
Tc such that at the critical temperature this mode goes over into the usual diffusive mode 
of the normal fluid. 

Our goal in the following is to establish the presence of the second sound pole in the 
holographic superconductor model [5]. In section two we introduce the model and describe 
its properties. We compute the condensate as a function of temperature. This is basically 
a review of the results of [5] except the fact that we choose to work in the grand canonical 

^Note that in [22] the mode has been called second sound whereas in [23] it was argued that it should 
be rather called fourth sound. There it has been argued that it is the fourth sound that survives the probe 
limit. In any case working directly in the probe limit we only have one sound mode. For convenience we 
will refer to it as second sound. Disentangling first, second or fourth sound would need to take into account 
the backreaction which is beyond the scope of this paper. 
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ensemble where we hold fixed the value of the chemical potential instead of the value of 
the charge density. 

In section three we compute the quasinormal modes of the complex scalar field in 
the unbroken phase. As expected we find that at the critical temperature a quasinormal 
frequency crosses over into the upper half of the complex frequency plane. Since a pole 
in the upper half plane is interpreted as an instability this is an indication that the scalar 
field condenses. The quasinormal modes of gauge fields in the four dimensional AdS black 
hole have been studied before [34]. In particular the longitudinal gauge field channel shows 
a diffusion pole with diffusion constant D = 3/(47rT). The holographic Green functions 
for gauge fields are often calculated in a formalism that employs gauge invariant variables, 

1. e. the electric field strengths. For reasons explained in section four we prefer however 
to work directly with the gauge fields. The longitudinal components obey two coupled 
differential equations subject to a constraint and we show in full generality how to compute 
the holographic Green function in such a situation. The quasinormal mode condition boils 
down to setting a determinant of field values at the boundary of AdS to zero. We find that 
this determinant is proportional to the electric field strength exemplifying that the poles 
of holographic Green functions are gauge invariant as expected on general grounds. 

Section four is the core of our paper. Here we study the lowest modes of the quasinormal 
mode spectrum in the superfluid phase. We find that the longitudinal gauge fields at non 
zero momentum couple to real and imaginary part of the scalar field fluctuations, building 
up a system of four coupled differential equations subject to one constraint. We solve this 
system numerically and compute the quasinormal modes from our determinant condition. 
We find hydrodynamic modes with approximately linear real part of the dispersion relation. 
We compute the speed of second sound from it and find our results to be in good numerical 
agreement with what was found in [22] from thermodynamic considerations. The second 
sound pole has however also an imaginary part and we can fit the dispersion relation (for 
small momenta) to a; = -izVgk — iVsk"^ which allows us to read off the attenuation constant 
Vg of second sound. We also find a purely imaginary mode with a dispersion relation of the 
form oj = —ij — iDk'^ with D the diffusion constant. It is a sort of gapped diffusion mode. 
The gap 7 goes to zero for T ^ Tc- A simple two fluid model suggests that there is still a 
normal fluid component and in it charges should diffuse in the usual way right below Tc. 
Therefore we expect a diffusive pole to show up in the two-point function. The diffusive 
behavior is modified however at long wavelength by the presence of the gap 7. This mode 
is therefore not really a hydrodynamic mode. 

We close this work with section five where we summarize and discuss our results. 

2. The Model 

As in [5] we consider a four dimensional planar AdS black hole with line element 
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the blackening factor is /(r) = — This metric has a horizon at rn = M^^^L"^^^, the 
Hawking temperature is T = In the fohowing we will rescale coordinates according 

to 

/ rHP \ 

L^THt 



(r\ 

t 

X 



\ L'^/ru y ) 



In the new dimensionless coordinates the metric takes the form (2J) with M 
the overall AdS scale I? . 

We take an abelian gauge model with a massive charged scalar field 



(2.2) 



1 times 
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and a tachyonic mass = —2/L? above the Breitenlohner-Priedmann bound. As in [5] 
we ignore the backreaction of these fields onto the metric. We seek solutions for which the 
time component of the gauge field vanishes at the horizon and takes a non-zero value p. 
on the boundary. This value can be interpreted as the chemical potential. The boundary 
condition on the horizon is usually justified by demanding that the gauge field has finite 
norm there. Here this can be seen as follows: as in [22] we can chose a gauge that removes 
the phase of the scalar field ^' from the equations of motion. In this gauge the scalar field 
current becomes = il^'A^. Therefore the value of is directly related to a physical 
quantity, the current, and it is a well defined physical condition to demand the current to 
have finite norm at the horizon. This is achieved by taking the scalar field to be regular 
and the gauge field to vanish at the horizon. 

In addition to the gauge field the scalar field might be non trivial as well. In fact for 
high chemical potential (low temperature) the scalar field needs to be switched on in order 
to have a stable solution. In our study of the quasinormal modes in the next section we 
will indeed see that a quasinormal mode crosses into the upper half plane at the critical 
temperature. We denote the temporal component of the gauge field in the dimensionless 
coordinates by ^. The field equations for the background fields are 
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^' = 0, 



-^> = 0. 



(2.4) 
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The equations can be solved numerically by integration from the horizon out to the bound- 
ary. As we just argued for the current to have finite norm at the horizon we have to 
chose $(1) = and demand the scalar field to be regular at the horizon. These conditions 
leave two integration constants undetermined. The behavior of the fields at the conformal 
boundary is 
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The value of the mass of the scalar field chosen allows to define two different theories due 
to the fact that both terms in the expansion above are normalizable in AdS. The canonical 
choice of what one considers to be the normalizable mode gives a theory in which -01 is 
interpreted as a coupling and '02 as expectation value of an operator with mass dimension 
two. On the other hand one might consider 02 as the coupling and ipi as the expectation 
value of an operator of dimension one. 

All our numerical calculations are done using the dimensionless coordinates. In order 
to relate the asymptotic values (^), ( ^.71) to the physical quantities we note that 
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(2.8) 
(2.9) 
(2.10) 
(2.11) 



where /i is the chemical potential, n the charge density and (Oi) are the vacuum expectation 
values of the operators sourced by the scalar field. Prom now on we will set L = 1 and 
work in the grand canonical ensemble by fixing // = 1. Different values for fl can now be 
interpreted as varying the temperature T. For high temperatures the scalar field is trivial 
and the gauge field equation is solved hy ^ = ft — ^. Spontaneous symmetry breaking 
means that an operator has a non trivial expectation value even when no source for the 
operator is switched on. We therefore look for nontrivial solutions of the scalar field with 
either = or -02 = 0. Numerically we find that a non-trivial scalar field is switched on 
at fl = 1.1204 corresponding to a critical temperature = 0.213^ for the operator Oi and 
at /2 = 4.0637 corresponding to a critical temperature of Tc = 0.0587^ for the operator O2. 
We chose to plot the squares of the condensates as a function of reduced temperature. It 
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Figure 1: The condensates as function of the temperature in the two possible theories. 



makes the linear behavior for temperatures just below the critical one manifest. 



(2.12) 



-6- 



Since we want to compute the holographic two point functions we will have to expand 
the action to second order in field fluctuations around the background. We divide the flelds 
into background plus fluctuations in the following way 



^ = V'(p) + (t{p, t, x) + iriip, t, x) , (2.13) 
A^, = Am{p) + a^{p, t, x) . (2.14) 



The gauge transformations act only on the fluctuations 

Saf, = df,X, (2.15) 
5a = -Xr], (2.16) 
6r] = Xa + Xij . (2.17) 

The action expanded out to second order is S = S^^^ + S^^^ + S^"^^ 
-5^°^ = / ^(-^V'''^ + ;^^'-(^mV')(5'^^)-^m-^^V'') , 

-5^ = J V=^(^-ijr^,/M- + lV-25^V'5'^a-2av,^''r?-2^V + 2^M5'^r?-2^^aV) , 

-2A''a^,il)a - 2df,aA''v + 2A^,adix'q + 2d„ria''ij - a^V'^ j . (2.18) 
Up to the equations of motion these can be written as boundary terms 
Sb^ = I V^{9'':Fp^a^-2gfPdp^a + 2APi;rj) , 

JB 

S^B =-J^ V^a"' Qs'Vp^aA + vdpV + (rdpa + a^V^) ■ (2-19) 

Note that S^^^ is not trivial since according to the holographic dictionary it has to encode 
the non-vanishing expectation values of the field theory operators. 

3. Quasinormal Frequencies in the Unbroken Phase 

We assume now that the scalar field ^ has vanishing background value and take the field 
to depend on t,x, p. Frequency uj and momentum k in the dimensionless coordinates are 
related to the physical ones ujp^, kp^ as 

u = ^ , k = f^. (3.1) 
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The equations of motion in the unbroken phase are 
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(3.2) 



The equation of motion for the complex conjugate scalar ^ can be obtained by changing 
the sign of the gauge field background $ in (p.2|). 



3.1 Green Functions 

In order to calculate the quasinormal frequencies we impose ingoing boundary conditions 
at the horizon. Since the coefficients of the differential equations ( |3.2D are known analyt- 
ically and are such that they are of Fuchsian type we can use the Frobenius method to 
approximate the solutions at the horizon and at the boundary by series expansions. As is 
well-known the holographic Green functions are proportional to the ratio of the connection 
coefficients. More precisely we demand 



H 



(p-1)— + 



on the horizon and write the local solution at the AdS boundary as 



B 



A B ^ 
- + ^ + 
P P^ 



P' 



(3.3) 



(3.4) 



In the theory with the dimension two operator we take A as the coefficient of the non- 
normalizable mode and B as the coefficient of the normalizable mode. Writing the local 
solution on the horizon as a linear combination of normalizable and non-normalizable modes 
on the boundary fixes the connection coefficient A and B. We have written the boundary 
action as a functional of real and imaginary part of the scalar field. We will rewrite the 
boundary action now in terms of the complex scalar ^ and its conjugate. In addition we 
introduce a local boundary counterterm to regularize the action. 



oB 



1 



(3.5) 



p=A 



This allows to compute the Green functions Cqq^q) = {0{—q)0{q)) and Gooil) = 
{0{—q)0{q)) fulfilling Goo{~Q) = GQQ{q), with q being the four momentum {uj,k). We 
write '${q,p) = '^o{Q)fq{p), where we interpret ^o('i') as the source that inserts the op- 
erator 0{q) in the dual field theory. We introduce a cutoff and normalize the profile 
function fg to 1/A at the cutoff. In terms of an arbitrary solution the normalized one is 
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fk{p) = where ^k{p) has the boundary expansion (|3.4| ). The boundary 

action is now 

= y^o(-'?)-^*^(A)^o(g) + ^o(-'?)-^*^(A)^o(g). (3.6) 

According to the holographic dictionary the renormahzed retarded Green functions are 
given by the hmit hniA^oo —2J^{A): 



G02O2 — ~^ ' ^0202 — ]| • (^•'^) 



We denote the connection coefficients for the complex conjugate scalar as B^. Accord- 
ing to ( p.2| ) they can be obtained by switching the sign of the chemical potential in the 
expressions for A,B. 

The theory with the operator of dimension one can be obtained through a Legendre 
transform. We note that the expectation value of the dimension two operator is (O2) = 
— p^{p"i>{p))'\^^^, whereas the source is given by A^'(A). We therefore add the following 
terms to the boundary action 

S^^S^ + J d'k [p\^{p^{p))' + ^{p^{p))'] 1^^^ . (3.8) 

Now we can evaluate the Green functions as before, the only difference being the normal- 
ization of the profile function fk{p) = ^fc(p)/[A-^''(A) -|- ^'(A)]. This normalization takes 
care that the term of order l/p^ in the boundary expansion couples with unit strength 
to the source ^^{q). We find finally for the Green functions of the Legendre transformed 
theory 

'^OiOi = ^ ' ^OiOi = ^ 5 (3-9) 

as expected. The quasinormal modes in the scalar sector are given by the zeroes of the 
connection coefficients A and A in the theory with operator of dimension two and by 
the zeroes of B and B in the theory with the dimension one operator. In terms of the 
unnormalized solutions to the field equations we can write the Green functions as 



Go^o2 = - [^'-^,+^] , (3.10) 



^0202 = - ^ 

^0,0, - A(AvI/;(A) + vI/,(A)) ■ ^-^-^^^ 
3.2 Quasinormal Modes from Determinants 

Before presenting the results for the quasinormal modes of the scalar field we would like to 
outline a method of how to calculate the holographic Green functions for the gauge fields 
without using gauge invariant variables such as the electric field strength E = —i(kat+ioax)- 



^Note that the infalhng boundary condition for the conjugate scalar is ^' ~ (p — 1) 
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The complicated structure of the gauge symmetry in the broken phase makes it rather 
difficuh to express the boundary action in terms of gauge invariant field combinations. As 
a warm up for the problem of how to calculate the holographic Green functions in this 
situation we will consider how we can calculate them in the unbroken phase directly in 
terms of the gauge fields. We necessarily have to solve a system of coupled differential 
equations whose solutions are restricted by a constraint. 

The correct boundary conditions for the gauge fields on the horizon are 



The two coefficients and are not independent but related by the constraint. At this 
point we have fixed the incoming wave boundary conditions and there seems to be now a 
unique solution to the field equations. We would expect however two linearly independent 
solutions with infalling boundary conditions on the horizon. The constraint reduces this 
to only one solution. The problem is now that in order to compute the Green function 
for the charge density and longitudinal current component separately we need solutions 
that asymptote to {at,ax) = (1,0) and {at,ax) = (0,1) respectively. This is of course not 
possible with only one available solution at the horizon. Because of the gauge symmetry 
the gauge field system (|3.2D allows for an algebraic solution 



with A' = 0, i.e. A being independent of p. This is of course nothing but a gauge transfor- 
mation of the trivial solution. Remember that even after fixing the radial gauge ap = 0, 
gauge transformations with gauge parameters independent of p are still possible. These 
gauge transformations appear as algebraic solutions to the field equations. We also stress 
that the infalling boundary conditions really have to be imposed only on physical fields, i.e. 
the electric field strength. Having therefore an arbitrary non trivial gauge field solution 
corresponding to an electric field with infalling boundary conditions we can add to it the 
gauge mode ( ^.14|) . 

We can use this to construct a basis of solutions that allows the calculation of the 
holographic Green functions. Let us now assume that there is a solution that takes the 
values {at, ax) = (1,0) at the boundary. We will call this solution from now on a- for 
i G t,x. Analogously we define the solution af. According to the holographic dictionary 
the solution a* couples to the boundary value limp^oo «*(<?, /o) = -^tiQ), i-e. the source of 
the field theory operator Jt (the time component of the conserved current J^). In parallel 
af couples to the boundary value Ax{q)- A generic solution of the gauge field equations 
can now be written in terms of the boundary fields as 




(3.12) 
(3.13) 



at = —ujX 
at = kX, 



(3.14) 
(3.15) 



ai{q,p) = Ax{q)af{p) + At{q)al{p) . 



(3.16) 
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Using this expansion the boundary action can be written as 

Sb = \j^M-Q) (^p'^o:l{-Q,p)-^ioiiiQ,p)) - f{pWx{-q,p)^iaiiq,p))^ Aj{q) , 

(3.17) 

where again we have introduced a cutoff at p = A. From this it follows that the holographic 
Green functions are given by 

2r^{p) = p^al{-k,p)^^{al{k,p))-f{p)al{-k,p)^^{aU^^ (3.18) 

in the limit 

G'^ = lim -2r^{A). (3.19) 

A— >oo 

Notice also that ^[.^*^'(p) - ^*^'{p)] = by the field equations. 

Although there are no terms of the form a[ax in the boundary action this formalism 
gives automatically expressions for the mixed Green functions Gtx and Gxt^- But we still 
have to construct the solutions a^j{p)- This can be done in the following way: suppose we 
have an arbitrary solution {at{p),ax{p)) obeying the infalling boundary conditions ( 3.12| ). 
We can add to this now an appropriate gauge mode, such that at the cutoff the solution 
takes the form (1,0) or (0, 1) in terms of at{A),ax{A). This is easily achieved by solving 




at (A) a^(A) 
—LuX kX 




(3.20) 



The linear combinations formed with the coefficients c*- give now new solutions (aj(/9), a^(yo)) 
and {af{p),ax{p))) obeying the correct boundary conditions on the AdS boundary. Using 
the general expression for the Green function ( 3.1S| ) we get explicitly in terms of solutions 
obeying the infalling boundary conditions 

G,, = lim A2--MA)^, (3.21) 

G,, = hm A^-—^^^ML—-, (3.22) 

Gxt = - lim A \ , (3.23) 

Gx. = - lim A \ ..^"^^"^^ . (3.24) 
A-.00 kat{A) + ujaxiA) ^ ' 

Note that the denominator for all is given by kat{A) +ujax{A) which is up to an irrelevant 
constant nothing but the gauge invariant electric field Ex- Therefore we see immediately 
that the poles of these Green function are gauge invariant and coincide of course with the 
poles of the Green function in the gauge invariant formalism where G (x Indeed using 



the constraint on the boundary we find the well known expressions [35] 



k"^ .oEl _ kio .oEi. 

tx 



^tt = 172 2 1™ ^^-^ ' = T2 2 J™ ^^-ET 

— A^oo Ex fc — W A->oo Ex 



,2 



Gx. = lim A^§^ . (3.25) 

fc^ — A^oo Ex 
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On general grounds one expects indeed that the poles of the holographic Green functions 
for gauge fields are gauge independent. 

If we are interested only in the location of quasinormal frequencies we do not even have 
to construct the holographic Green functions explicitly. From the linear system in ( 3.2C| ) 



we infer that the quasinormal frequencies coincide with the zeroes of the determinant of the 
field values at the boundary. Indeed vanishing determinant means that there is a nontrivial 
zero mode solution to ( |3.20D such that the boundary values of the fields are (0, 0) which in 
turn means that the coefficient in the solution of the non-normalizable mode vanishes. The 
determinant is X{kat + loUx) and again given by the electric field strength. In fact these 
remarks apply to systems of coupled differential equations in AdS black hole metric in 
general: the quasinormal frequencies corresponding to the poles of the holographic Green 
functions are the zeroes of the determinant of the field values on the boundary for a maximal 
set of linearly independent solutions obeying infalling boundary conditions on the horizon. 
The fact that the differential equations are coupled is the holographic manifestation of 
mixing of operators under the RG flow. Therefore one has to specify at which scale one is 
defining the operators. The scheme outlined above is dual to define the operators at the 
cutoff A. 

3.3 Hydrodynamic and higher QNMs 

We have numerically computed the quasinormal frequencies for the fluctuations satisfying 



the equations of motion (|3.2| ) for both the O2 and the Oi theories. The quasinormal modes 
of the scalar field correspond to zeroes of A in the theory of dimension two operator and 
to zeroes of B in the dimension one operator theory, where A and B are the connection 
coefficients of the boundary solution (p.4|). Results for the lowest three poles of the scalar 
field at zero momentum are shown in Figure ^. 

The poles with positive real part correspond to the quasinormal modes of the complex 
scalar ^, while those with negative real part are the quasinormal modes of ^, obtained 
by changing the overall sign of the gauge field background As the temperature is 
decreased, the poles get closer to the real axis, until at the critical temperature Tc the 
lowest mode crosses into the upper half of the complex frequency plane. It happens at 
Tc = 4.0637 for the theory of dimension two operator and at Tc = 1.1204 in the case 
of the dimension one operator theory. For T < Tc the mode would become tachyonic, 
i. e. unstable. This instability indicates that the scalar field condenses and the system 
undergoes a phase transition at T = Tc. At the critical point, the lowest scalar quasinormal 
mode is a hydrodynamic mode in the sense that it is massless, limk-^QUj{k) = 0. This mode 
is identified with the Goldstone boson that appears after the spontaneous breaking of the 
global U{1) symmetry and in the next section we will see that it evolves into the second 
sound mode characteristic of superfluid models. 

The quasinormal modes correspond to simple poles of the retarded Green function, so 
close to the nth pole the Green function can be approximated by G{uj,k,T) ~ l^l^'^f^^) • 
Knowing the connection coefficients we can compute the Green functions and therefore the 
residue for each quasinormal mode as explained in [36,37]. For the lowest quasinormal 
mode at A; = and at the critical temperature, the residue takes the value R2{Tc) = 
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Figure 2: Lowest scalar quasinormal frequencies as a function of the temperature and at momen- 
tum k = 0, from T/T^ = 00 to T/Tc = 0.81 in the O2 theory (right) and to T/T^ = 0.56 in the Oi 
theory (left). The dots correspond to the critical point T/Tc — 1 where the phase transition takes 
place. Red, blue and green correspond to first, second and third mode respectively. 



-2.545 + 0.825i in the O2 theory and i?i(Tc) = 0.686 -0.348i in the Oi theory. In general, 
one expects the residues of hydrodynamic modes that correspond to conserved quantities 
of the system to vanish in the limit of zero momentum, since its susceptibility remains 
constant. Consider for instance the diffusion mode associated to conserved density. The 
susceptibility is defined through the two point correlation function as 

iak'^ a 

X = lim (pp) = lim — — r = — , (3.26) 

where D is the diffusion constant and a is the conductivity. The residue, iak'^, vanishes and 
one recovers the well-known Einstein relation a = Dx- However, for hydrodynamic modes 
appearing at second order phase transitions the order parameter susceptibility should di- 
verge at the critical point. This order parameter susceptibility is given in our case by the 
correlator of the boundary operator sourced by the scalar field. At the critical temperature 
it is 

Xo.o. = lim (O.O.) = hm - 00 (3.27) 

* * k,LL}^0 k,ui-*0 UJ — LOH[k, Ic) 

since a;//(0, Tc) = while the residue remains finite. This result allows us to identify the 
lowest scalar quasinormal mode in the unbroken phase with the Goldstone boson appearing 
at the critical point. 

In the model under consideration one can also compute the gauge field fluctuations 
in the normal phase. Nevertheless, as the model does not include the backreaction of 
the metric, the computation is not sensitive to temperature anymore. This can be seen 
from the equations of motion (|3.2| ) of the gauge fluctuations, that do not depend on the 
background solutions thus are independent of the temperature. Hence we recover the 
results for the quasinormal modes of vector field perturbations in the AdS4 black hole 
background computed by [34]. For our purposes the most important fact is the presence 



-13- 



of a hydrodynamic mode corresponding to diffusion. For small momenta that mode has 
dispersion relation uj = —iDk'^ with D = \ (which is -D = 3/(47rr) in physical units). In 
order to study the behavior of the diffusion pole in the unbroken phase as a function of the 
temperature one has to consider the backreacted model described in [38]. 



4. Quasinormal Frequencies in the Broken Phase 

In this section we will apply our determinant method for finding quasinormal modes of a 
coupled system of field equations. With this technique we follow the model analyzed in the 
previous section into its broken phase. Figure |3| schematically summarizes our analysis. 
The two formerly separate sets of scalar (grey dots) and longitudinal vector poles (black 
dots) present in the unbroken phase where the scalars and vectors decouple, are now unified 
into one inseparable pole structure in the coupled system. This is in analogy to a coupled 
system of two harmonic oscillators in which it makes no sense to ask for the eigenfrequencies 
of the single oscillators. One could of course try to diagonalize the system of differential 
equations, in our case however this looks rather complicated and we prefer to work directly 
with the coupled system and with the gauge fields instead of gauge invariant variables. The 
lowest modes (see figure ^) are two hydrodynamic second sound modes originating from 
the two lowest scalar quasinormal modes in the unbroken phase. In addition we will find 
a non-hydrodynamic pseudodiffusion mode staying on the imaginary axis in the range of 
momenta we consider. This mode can be thought of as the prolongation of the diffusion 
mode into the unbroken phase. 

4.1 Application of the determinant method 

The equations of motion in the broken phase couple the scalar fluctuations r/, a to the 
longitudinal vector components at, ax 

= fv + /H ] V + [ -r + + 9P/ 



(4.1) 

„ = V.,(;',_i),.,(^,_,___),,^<., + _*,. ,4.2, 

= /a/' + —a/ - f ^ + 2V^2^ _ _ 2iu:^ rj-AiPcPa, (4.3) 

P VP / 

= fa/ + faj + (y - 2^^') ax- + y «* + 2ik^lj r/ . (4.4) 

This system of four coupled equations is subject to the constraint 

^a/ + \aj = 2i{iP'r]-ij r/') , (4.5) 

where the left hand side known from the unbroken phase is amended by the condensate 
terms on the right. Note that the real part a of the scalar fluctuation is not involved in 
the constraint. 
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Figure 3: Schematic plot of the poles in the coupled system, i.e. in the broken phase at small 
finite momentum right below Tc. These poles are present in each retarded correlation function 
for the coupled fields rj, a, At, A^, while their residues might vanish for specific fields. Close to 
the origin we find the (pseudo)diffusion mode and two hydrodynamic second sound modes. In 
addition two sets of higher (non-hydrodynamic) quasinormal modes are shown. In the unbroken 
phase these poles originate from the scalar (grey dots). We also expect a tower of purely imaginary 
poles stemming from the longitudinal vector channel (black dots). The grey area indicates where 
our present numerical methods break down. In this paper we will be concerned primarily with the 
hydrodynamic modes and will touch upon the higher quasinormal modes only briefly. 



The constraint can be interpreted as the Ward identity of current conservation in the 

presence of the condensate. We expand the gauge fields near the boundary and note that 
Oq = {n) / and = —{jx)/p'^, where jx is the x-component of the current. Expanding 
also the r.h.s. and comparing the leading orders in p we find 

d^{f) = 2{Oi)r^i (4.6) 

where rf^ is the source for the insertion of the imaginary part of the operator Oj, i.e. the 
Goldstone field in the dual field theory. This equation is to be understood as the local 
Ward identity encoding current conservation in the presence of the condensate (Oj). It 
follows for example that the two point function of the divergence of the current with the 
operator O^^^ is zero only up to a contact term {d^j'^(x)ol^\y)) = {Oi)S{x — y). 

The gauge field component ay being transverse to the momentum decouples from the 
above system and assumes the form 

= / a/ + fay' +(^^-^-2^P^^ ay. (4.7) 

Since we do not expect any hydrodynamic modes in the transverse vector channel we will 
not study this equation further. 
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Applying the indicial procedure to the system ( [4 .11) to find the exponents for the 
singular and the coefficients for the regular parts of the fields, we obtain the following 
behavior at the horizon 



V={p-lf (r/W+r?«(p-l) + ...) , 
{p-lf (a(°)+a«(p-l) + ...) , 

(p-l)C (ajna«(p-l) + ...' 



a 



at 



(4.8) 
(4.9) 
(4.10) 
(4.11) 



with the exponent C = —ioo/'i obeying the incoming wave boundary condition. 

Due to the constraint we can choose only three of the four parameters at the horizon. 
Using the constraint and without loss of generality we can eliminate the time component 
and parametrize the solutions by {r]'^^\ We choose three linearly independent 

combinations /, //, ///. A fourth solution can be found from the gauge transformations 



7/ 



IV 



IV 



0. 



jv 



JV 



-Xk. 



(4.12) 



with A being an arbitrary constant with respect to p. It is not an algebraic solution to the 
equations of motion since rj has non trivial dependence on the bulk variable p. The gauge 
solution solves the equations (^]|) not exactly but only up to terms proportional to the 



background equations (2.4). 

Our goal is to find the poles in the retarded correlation functions of the four fields 
appearing in the coupled system of equations of motion (|4.1| ). A convenient way of imposing 
the appropriate boundary conditions is given by redefining the scalar fields as 



flip) = PVip) 



d{p) = pa{p) 



(4.13) 



Then the most general solution for each field (pi G {77, a, at, ax} including gauge degrees 
of freedom can be written 



'Pi 



(4.14) 



In the theory with the dimension two operator the sources for the various gauge in- 
variant operators are given by ipi{K). We are interested in the quasinormal modes of the 
system ( [4.1[ ) and as we have argued in the previous section these are the special values of the 
frequency where the determinant spanned by the values ' ' vanishes. Expanding 
this determinant we get 



^det 
A 



pj' vj" pj"" 



d' 
dJ 



dt" 



III 



d^ 
d^ 



(4.15) 



i^rf det 



pj fJ^ pJ^^ 



d' 

X dJ 



+ uj det 



\ dx 



dJ 



III\ 
) 



+ A;det 



iii\ 
III 



JI 
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which needs to be evaluated at the cutoff p = A. The first term in ( 4.15 ) vanishes at the 
cutoff since (p^ = Kifj = is just the condition that the operator O2 is not sourced by the 
background. 

We first find three hnearly independent numerical solutions and then solve condi- 
tion (4.15) numerically. Explicit checks confirm that the choice of a solution basis ^p^'^^'^^^'^^ 



is completely arbitrary and does not change the results. Note also that in our present case 
all the remaining determinants can not be factorized. But if the momentum is set to zero, 
the only remaining term is the one with u) and the determinant factorizes into a scalar part 
and a vector part since the system of equations decouples. 

In the theory with the dimension one operator the sources are given by —A^fj and 
— A^CT for the scalar fields. Let us call <pi = —p^fj' and (^2 = —p^cr' in this case. The 
determinant has therefore the same form and again the first term vanishes due to the 
absence of sources for Oi in the background solution. The quasinormal modes can again 
be found by integrating three arbitrary solutions with infalling boundary conditions from 



the horizon to the cutoff and finding numerically the zeroes of the determinant ( 4.15 ). 
4.2 Hydrodynamic and Goldstone modes 

Sound mode The scalar modes originally destabilizing the unbroken phase turn into 
Goldstone modes at Tc instead of becoming tachyonic. Below Tc they evolve into the 
two second sound modes. Figure ^ shows their movement when momentum is changed at 
different temperatures. Note that we focus on the positive real frequency axis because of the 
mirror symmetry sketched in figure ^. From the dispersion relation at small frequencies 
and long wavelengths we extract the speed of second sound Vg and the second sound 
attenuation F^ using the hydrodynamic equation 

uj = Vsk-iTsk'^ . (4.16) 



It turns out that the hydrodynamic regime, i.e. that range of momenta in which the 
dispersion relation is well approximated by (4.16), is very narrow for temperatures just 



below the critical one since the speed of sound vanishes at Tc- Fits to the hydrodynamic 
form at a high temperature T ~ O.OOOOTc are plotted in figure [s] for the 02-theory, the 
results for the Oi theory are qualitatively similar. 

The speed of second sound is shown in figure ^. We have a good numerical agreement 
with the thermodynamic value of the second sound velocity given in [22]. This nicely 
confirms validity of our method. In particular, within our numerical precision we find that 
the value of the square of the speed of sound tends to tig 1/3 in the Oi theory and 
to v1 ~ 1/2 in the O2 theory'^ . However, the numerics becomes rather unstable for low 
temperatures, especially for the Oi theory. Near but below the critical temperature we 



^In [23] it was argued that conformal symmetry implies v1 — 1/2 at zero temperature. Due to the 
divergence in the order parameter for the Oi theory conformal symmetry could be broken and allow thus 
for a different value. 
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Figure 4: Movement of the positive frequency sound pole away from uj — with increasing 
spatial momentum. Distinct curves correspond to temperatures below the phase transition T/Tc = 
0.999 (black), 0.97 (pink), 0.91 (red), 0.71 (green), 0.52 (blue), 0.26 (light blue). Dots on one curve 
are separated by Afc — 0.05. All curves start at fc = 0.05 and end at fc = 1.00. 




Figure 5: Fits of the real and imaginary part of the hydrodynamic modes in the broken phase to 
the lowest order approximation uj — Vgk — iTsk^- The left figure shows the real part and the right 
one the imaginary part. The thick lines are the numerical results and the thin lines are the linear 
and quadratic fits. The fit is done for a temperature just below the critical one where the range of 
the approximation is rather small. 

find 

01- Theory, (4.17) 

02- Theory. (4.18) 

Moreover, as a benefit of our effort considering the fluctuations, we are also able to 
extract non-thermodynamic quantities in this channel. Specifically we examine the atten- 
uation of the second sound mode as shown in figure |^. The curve shows how attenuation 
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Figure 6: The plots show the squares of the speed of sound as extracted from the location of the 
lowest quasinormal mode in the broken phase. The left figure is for the Oi theory and the right one 
for the O2 theory. We also indicate the linear behavior close to Tc. As can be seen the numerics 
for the Oi theory becomes somewhat unstable for low temperatures. 

smoothly asymptotes to zero as the superfluid becomes more and more ideal at low tem- 
peratures. This effect is however much stronger in the O2 theory. Near Tc the attenuation 
is growing. Within our numerical precision it seems however that the attenuation constant 
is taking a finite value at the critical temperature. Numerically we find 

r, = l.STTc at r = 0.999irc Oi- Theory, (4.19) 
= 1.48rc at r = 0.9998Tc O2- Theory. (4.20) 

A similar behavior has been observed in [39], where the attenuation of the normal sound 
mode asymptotes to a finite value near a phase transition. 




I — 0.0 1 — ' — ' — ' — < — ^ I — ' — < — ' — ' — ' — ' — 
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Figure 7: The plots show the attenuation constants of second sound as extracted from the location 
of the lowest quasinormal mode in the broken phase. The left figure is for the Oi theory and the 
right one for the O2 theory. Again it can be noticed that the Oi theory is numerically more 
challenging at low temperatures. 

(Pseudo) diffusion mode The vector diffusion mode from the unbroken phase turns 
into a (pseudo) diffusion mode below Tc^. For not too low temperatures and not too large 

*An analytical result obtained for second sound in a non-abelian model [24] also shows the appearance 
of a pseudo diffusion mode with a gap that vanishes as the condensate goes to zero. 
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momenta the dispersion relation for this mode is weh approximated by 

uj = -iDk^ - i-f{T) , 



(4.21) 



with a gap 7 E M in imaginary frequency direction. Thus the pole is shifted from its 
unbroken phase position such that it does not approach zero at vanishing momentum any 
more, i.e. it is not anymore a hydrodynamic mode. 




0.85 0.90 0.95 1.00 T,. "0.7O 0.75 0.80 0.85 0.90 0.95 1.00 7;. 

Figure 8: The plots show the gap of the pseudo diffusion mode as a function of reduced tempera- 
ture. On the left the Oi theory and on the right the O2 theory. 



In figure ^ we have plotted the gap 7 as a function of the reduced temperature and we 
can see that it vanishes linearly near Tc. 




Figure 9: The plots show the dispersion relations for the pseudo diffusion mode at dilferent tem- 
peratures and the diffusion dispersion relation with D — 3/(47rT) (note that even in the unbroken 
phase the latter approximates the diffusive quasinormal mode only for small momenta). On the 
left we have the Oi theory at temperatures T = 0.999rc (black), T = O.QTTc (red) and T = O.OlTc 
(green). On the right the same for the O2 theory at temperatures T = 0.999Tc (black), T = 0.97Tc 
(red) and T = O.STTc (green). 

Figure ^ shows the dispersion relation for the diffusion pole at different temperatures. 
The offset at = is the gap size 7 depending linearly on T only near Tc- This implies 
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that the relation ( 4.21 ) asymptotes to the ordinary diffusion equation near the critical 
temperature. As expected the highest temperature curve T = 0.999Tc (black) matches 
the hydrodynamic approximation (thin line) very well at small momenta. That agreement 
becomes worse around k ~ 0.25. Also as the condensate grows below the behavior of 
this (pseudo) diffusion mode becomes less hydrodynamic. 



4.3 Higher quasinormal modes 

In addition to the hydrodynamic sound modes and the pseudo diffusion mode there are 
higher quasinormal modes. They are not the main focus of this paper and we have not 
studied them in detail. We have however traced the prolongation of the second and third 
quasinormal modes in the scalar sector from the unbroken phase into the broken phase. 

The former scalar modes evolve continuously into higher modes of the coupled system 
through the phase transition as seen from the two kinked dashed (unbroken phase) and solid 
(broken phase) lines in figure 10. The kink indicates that the poles move continuously but 
change direction at the critical temperature. We show only the plot for the theory with the 
dimension two operator, similar results hold for the Oi theory. At the critical temperature 
the locations of the quasinormal frequency calculated with the Frobenius method in the 
unbroken phase and by the method of finding the zeroes of the determinant spanned by 
the solutions match with impressively high precision. We might take this as a highly non 
trivial test of the accuracy of the numerical integration method. 

We expect that all quasinormal frequencies are shifted continuously in the complex 
frequency plane across the phase transition. This means that there are no jumps in any of 
the dispersion relations. There is simply an infinite set of poles corresponding to the degrees 
of freedom of the system which are continuously shifted when parameters are changed. 

Our numerical invesigations also reveal poles developing an increasing real part with 
decreasing temperature which most likely originate from the longitudinal vector modes in 
the unbroken phase. Note that we have not shown these modes for simplicity. 

The two higher modes shown in figure 10 move parallel to the real axis and to each 
other with decreasing temperature. ^ At low temperatures their real parts are almost the 
same. This is also true for the longitudinal vectors as far as we can tell within our numerical 
uncertainties. We suspect that this alignment of higher poles has to do with the appearence 
of the conductivity gap showing up at low temperatures [5]. 



5. Summary and Discussion 

We have studied the hydrodynamics of a holographic model of a superfluid. The model is 
an Abelian gauge field model with a charged scalar field in a four dimensional AdS black 
hole background. 

As is well known by now the holographic dictionary allows to interpret the quasinormal 
frequencies as the poles of the retarded Green functions of the dual field theory. By 

^Note that this behavior is very different from the behavior usually found in five-dimensional holographic 
setups when temperature is decreased [40-42]. 
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Figure 10: Movement of the higher poles in the complex frequency plane at vanishing momentum 
k = 0. Scalar modes 2 (blue) and 3 (green) in the unbroken phase (dashed) evolve continuously into 
the higher poles of the broken phase (solid). The right end points are evaluated for T/Tc = 0.25. 

calculating the low lying quasinormal frequencies numerically in the broken and unbroken 
phases we were able to identify the hydrodynamics. We found that at high temperatures, 
T > Tc there is only one hydrodynamic mode representing the diffusion of the conserved 
U{1) charge. As one lowers the temperature reaching the critical temperature Tc two of 
the quasinormal frequencies of the charged scalar field approach the origin of the complex 
frequency plane. Precisely at the critical temperature at the onset of the phase transition 
these modes become massless, giving rise to new hydrodynamic variables. We also have 
calculated the residue of these modes and found that it stays finite at T = Tc resulting in 
a divergence of the order parameter susceptibility, as expected. 

Below the critical temperature these modes stay massless and show a dispersion relation 
with a linear real part and a quadratic imaginary one, allowing an interpretation as the 
modes of second sound in the superfluid. On the other hand the diffusion mode starts to 
develop a gap and stops to be hydrodynamic. The counting is therefore one hydrodynamic 
mode at high temperatures, three at the critical temperature and two at low temperatures. 

In the low temperature phase we were able to calculate the speed of sound as well as 
its attenuation constant as a function of temperature. Also the gap in the pseudo diffusion 
mode has been determined. 

We have also been able to follow some of the higher quasinormal modes through the 
phase transition and found that they evolve continuously albeit non-smoothly with tem- 
perature in the complex frequency plane, showing a sharp kink at the critical temperature. 

On a technical side we have developed a method to determine the quasinormal frequen- 
cies and the holographic Green functions for systems of coupled differential equations and 
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without using gauge invariant variables. The quasinormal frequencies correspond simply to 
the zeroes of the determinant spanned by a maximal set of linearly independent solutions. 
We have furthermore seen that the poles of the Green functions stemming from bulk gauge 
fields are gauge invariant as expected. 

There are now several interesting questions that should be investigated in the future. 
The most obvious one is to apply the methods developed here to the model where the 
backreaction of the matter and gauge field onto the metric are properly taken into account. 
Due to the presence of the metric fluctuations the hydrodynamics of such a model is 
certainly much richer. In addition to the diffusion and second sound modes found here 
one expects also shear and sound modes stemming from bulk metric fluctuations. The 
corresponding system of differential equations promises to be rather involved. However 
there should be no principal obstacle to apply our methods also in these cases. Such an 
investigation is currently underway [43]. 

Another interesting direction of research should be to reinterpret the results obtained 
here and analogous ones for related models with different scalar mass and living in different 
dimensions from the point of view of dynamical critical phenomena [44]. We have seen 
already that the speed of sound scales with exponent one half, whereas the gap in the 
pseudo diffusion mode scales with exponent one. The situation for the sound attenuation 
is unfortunately less clear. As far as our numerics indicates the sound attenuation reaches 
a finite value at T = Tc. An extensive study of related models possibly with enhanced 
numerical efforts might give new insights here. 

Of course an effort should also be undertaken to extend the results at hand to models 
of p-wave superconductors [45]. For infinitesimal condensates and analytic study of the 
hydrodynamics has already been done in [24]. It might be of interest to supplement these 
analytical result with a numeric study that allows to go further away from the phase 
transition point deep into the broken phases. Another very interesting class of holographic 
p-wave superconductors are the ones realized on D7 brane embeddings [8,13,21]. Due to 
the presence of fundamental matter these should be especially interesting to study. 

We hope to come back to all or at least some of these and other questions in future 
publications. 
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